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Abstract 

We summarize the current status of our understanding of nuclear forces based upon chiral symmetry — 
an idea that was advocated by Gerry Brown already several decades ago. 

1 Introduction 

Forty years ago, Gerry Brown published the Comment "Isn't it Time to Calculate the Nucleon-Nucleon 
Force" [T, and more than 30 years ago, he wrote a book chapter entitled "Chiral Symmetry and the Nucleon- 
Nucleon Interaction" [2^. In fact, as early as 1968, he had published, together with two co-workers, a paper [3] 
on three-nucleon forces, where the consequences of chiral symmetry were fully exploited. Nevertheless, it 
should take a few more decades until the problem of the nucleon-nucleon (NN) interaction was really solved 
taking chiral symmetry consistently into account — and this fact proves the very advanced nature of Gerry's 
ideas concerning nuclear forces. It is the purpose of this contribution to summarize where we are today in 
our understanding of nuclear forces, based upon concepts that Gerry advocated already many decades ago. 

But let's first recall more of the interesting history of the idea of chiral symmetry. The modern un- 
derstanding is that this symmetry arises because the up and down quarks happen to have relatively small 
masses. However, chiral symmetry and its significance for low-energy hadron (pion) physics was discovered 
long before QCD. In 1960, based upon concepts proposed by Schwinger 14 , Gell-Mann and Levy [5] devel- 
oped the sigma model, which is a linear realization of chiral symmetryPJ One major problem researchers 
had been struggling with in the 1950's was that the pion-nucleon scattering length came out two orders of 
magnitude too large when the (renormalizable) pseudo-scalar (75) TriV interaction was used. This unrealistic 
prediction was due to very large contributions from virtual anti-nucleon states (the so-called "pair terms" 
or "Z-graphs"). Similar problems occurred in the 27r-exchange contribution to the NN interaction. In the 
sigma model, the large pair terms are cancelled by processes involving the (fictitious) cr boson. In this way, 
the linear sigma model demonstrates how imposing chiral invariance fixes the problem with low-energy tt-N 
scattering. However, the fictitious character of the a particle as well as the reliance on the perfect cancelation 

*Chapter 14 of: From Nuclei to Stars — Festschrift in Honor of Gerald E. Brown, edited by Sabine Lee (World Scientific, 
Singapore, 2011), p. 317-343. 
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of huge terms are uncomfortable features. In 1967, motivated by the current algebra approach to soft pion 
physics, Weinberg [7] worked out what has become known as the non-linear sigma model, which does not 
include a a anymore and has pions and nucleons interact via pseudo- vector (derivative, 757'^ coupling 
besides a new (non-linear) -kttNN term also involving a derivative ( "Weinberg- Tomozawa term" [8i j9j ) . The 
derivative (equivalent to momentum) guarantees that the interaction vanishes when the momentum goes to 
zero providing a natural explanation for the weakness of the interaction by soft pions which does not rely 
on the cancellation of large terms. Following suggestions by Schwinger, Weinberg [lOj developed, soon after, 
a general theory of non-linear realizations of chiral symmetry, which was further generalized in an elegant 
way by Callan, Coleman, Wess, and Zuimino [11 . 

However, ideas were needed for how to implement chiral symmetry consistently in the theory of pionic and 
nuclear interactions and how to deal with renormalization, since the derivative coupling is not renormalizable 
in the conventional sense. In his contribution to the 'Festschrift' in honor of Schwinger of 1979 [T2 j I13 j . 
Weinberg proposed to consider the most general possible Langrangian including all higher-derivative terms 
that are consistent with chiral symmetry (besides the other commonly assumed symmetry principles). For 
this theory to be manageable, one needs to assume some sort of perturbative expansion such that only a finite 
number of terms contribute at a given order. This expansion is provided by powers of small external momenta 
over the chiral symmetry breaking scale, A;^ ^ 1 GeV. The higher-derivative terms supply the counterterms 
that make possible an order-by-order renormalization, which is the appropriate renormalization procedure 
for an effective field theory. Weinberg's suggestions were soon picked up by Gasser, Leutwyler, and associates 
who worked out, to one loop, the cases of tttt [14^ and ttN scattering [TJ] with great success. 

But there was still the problem of the nuclear force which is more difficult, since nuclear interactions do 
not vanish in the chiral limit {q — > 0; mu/j_,mT^ — >■ 0) and require a non-perturbative treatment because of 
the existence of nuclear bound states. In a series of papers published around 1990 [THl HZl [TB], Weinberg 
picked up the nuclear force issue and suggested to calculate the NN potential perturbatively in the chiral 
expansion and then iterate it to all orders in a Schroedinger or Lippmann-Schwinger equation to obtain the 
nuclear amplitude. Here, the introduction of four-nucleon contact terms is crucial for renormalization. 

Following the Weinberg proposal, pioneering work was performed by Ordoiiez, Ray, and van Kolck |19|,l20j 
who applied time-ordered perturbation theory to construct a NN potential up to next-to-next-to-leading 
order (NNLO). The results were encouraging and nuclear EFT quickly developed into one of the most popular 
branches of modern nuclear physics. The Munich group used covariant perturbation theory and dimensional 
regularization to calculated the perturbative NN amplitude without [21| and with A(1232)-isobar degrees 
of freedom [22] at NNLO. Besides this, the Munich group worked out important loop contributions of 
higher order [23l [Ml [SSI [ZSl [27| . A relativistic approach was also taken by the Brazil group [28[ [29] . The 
Bochum-Jiilich group devised a method of unitarity transformations to eliminate the energy-dependence of 
time-ordered perturbation theory amplitudes and calculated the NN potentials up to NNLO [501 [31] ■ The 
Idaho group managed to construct a chiral NN potential at next-to-next-to-next-to-leading order (N^LO) 
and showed that only at this order can one achieve the precision necessary for reliable few-nucleon and 
nuclear structure calculations [32l [33l [34l [35] . Progress extended beyond the NN interaction, as nuclear 
many-body forces based upon chiral perturbation theory were also developed [HI [36l [ST] [38l [39] . 

During the past decade or so, chiral two-nucleon forces have been used in many microscopic calculations 
of nuclear reactions and structure [ini[M[ll[13[M[13[li[171[li[ini[5ni[SIl[Sl[Sl[51E^ and the 
combination of chiral two- and three-nucleon forces has been applied in few-nucleon reactions ^TJ [5^1 [53 [SS] 
[59, .60, 61] [621 [631 [64], structure of light- and medium-mass nuclei [65ll66, 67, 68 , and nuclear and neutron 
matter [69] [70] — with a great deal of success. The majority of nuclear structure calculations is nowadays 
based upon chiral forces. 

This article is organized as follows. In Section [2] we sketch the essential ideas of an EFT for low-energy 
QCD. Section [3] provides an overview on nuclear forces derived from chiral EFT. The two-nucleon force is 
then discussed in detail in Section [4] Many-body forces are the subject of Section [5] Finally, Section [6] 
contains our conclusions. 
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2 Low-energy QCD and effective field theory 



Quantum chromodynamics (QCD) is the theory of strong interactions. It deals with quarks, gluons and 
their interactions and is part of the Standard Model of Particle Physics. QCD is a non-Abelian gauge field 
theory with color SU (3) the underlying gauge group. The non-Abelian nature of the theory has dramatic 
consequences. While the interaction between colored objects is weak at short distances or high momentum 
transfer ("asymptotic freedom"); it is strong at long distances ( ~ 1 fm) or low energies, leading to the 
confinement of quarks into colorless objects, the hadrons. Consequently, QCD allows for a perturbative 
analysis at large energies, whereas it is highly non-perturbative in the low-energy regime. Nuclear physics 
resides at low energies and the force between nucleons is a residual color interaction similar to the van der 
Waals force between neutral molecules. Therefore, in terms of quarks and gluons, the nuclear force is a 
very complicated problem that, nevertheless, can be attacked with brute computing power on a discretized, 
Euclidean space-time lattice (known as lattice QCD). Advanced lattice QCD calculations [7T1[72] are under 
way and will yield improving results in the near future. However, since these calculations are very time- 
consuming and expensive, they can only be used to check a few representative key-issues. For everyday 
nuclear structure physics, a more efficient approach is needed. 

The efficient approach is an effective field theory. For the development of an EFT, it is crucial to identify 
a separation of scales. In the hadron spectrum, a large gap between the masses of the pions and the masses 
of the vector mesons, like p(770) and w(782), can clearly be identified. Thus, it is natural to assume that 
the pion mass sets the soft scale, Q ~ m^, and the rho mass the hard scale, A^^ ~ nip, also known as the 
chiral-symmetry breaking scale. This is suggestive of considering an expansion in terms of the soft scale 
over the hard scale, Q/A^. Concerning the relevant degrees of freedom, we noticed already that, for the 
ground state and the low-energy excitation spectrum of an atomic nucleus as well as for conventional nuclear 
reactions, quarks and gluons are ineffective degrees of freedom, while nucleons and pions are the appropriate 
ones. To make sure that this EFT is not just another phenomenology, it must have a firm link with QCD. 
The link is established by having the EFT observe all relevant symmetries of the underlying theory. This 
requirement is based upon a 'folk theorem' by Weinberg jl2| : 

If one writes down the most general possible Lagrangian, including all terms consistent with 
assumed symmetry principles, and then calculates matrix elements with this Lagrangian to any 
given order of perturbation theory, the result will simply be the most general possible S-matrix 
consistent with analyticity, perturbative unitarity, cluster decomposition, and the assumed sym- 
metry principles. 

Since the up and down quark masses are very small, one may assume (approximate) chiral symmetry 
which, however, is broken in two ways. It is explicitly broken, because the the up and down quark masses 
are not exactly zero. Moreover, we are also faced with so-called spontaneous symmetry breaking. 

A (continuous) symmetry is said to be spontaneously broken if a symmetry of the Lagrangian is not 
realized in the ground state of the system. There is evidence that the (approximate) chiral symmetry of the 
QCD Lagrangian is spontaneously broken — for dynamical reasons of nonperturbative origin which are not 
fully understood at this time. The most plausible evidence comes from the hadron spectrum. From chiral 
symmetry, one would naively expect the existence of degenerate hadron multiplets of opposite parity, i.e., 
for any hadron of positive parity one would expect a degenerate hadron state of negative parity and vice 
versa. However, these "parity doublets" are not observed in nature. For example, take the p-meson which is 
a vector meson of negative parity (J^ = 1^) and mass 776 MeV. There does exist a 1+ meson, the ai, but 
it has a mass of 1230 MeV and, thus, cannot be perceived as degenerate with the p. On the other hand, the 
p meson comes in three charge states (equivalent to three isospin states), the p^ and the p^, with masses 
that differ by at most a few MeV. In summary, in the hadron spectrum, SU{2)v (isospin symmetry) is well 
observed, while axial symmetry is broken: SU{2)l x SU{2)ii is broken down to SU{2)v- 

A spontaneously broken global symmetry implies the existence of (massless) Goldstone bosons [73l [74] . 
The Goldstone bosons are identified with the isospin triplet of the (pseudoscalar) pions, which explains why 
pions are so light. The pion masses are not exactly zero because the up and down quark masses are not 
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exactly zero either. Thus, pions are a truly remarkable species: they reflect spontaneous as well as explicit 
symmetry breaking. Goldstone bosons interact weakly at low energy. They are degenerate with the vacuum 
and, therefore, interactions between them must vanish at zero momentum and in the chiral limit (771^^ 0). 

The next step is to build the most general Lagrangian consistent with the (broken) symmetries discussed 
above. An elegant formalism for the construction of such Lagrangians was developed by Callan, Coleman, 
Wess, and Zumino (CCWZ) [TT] who worked out the group-theoretical foundations of non-linear realizations 
of chiral symmetry. It is characteristic for these non-linear realizations that, whenever functions of the 
Goldstone bosons appear in the Langrangian, they are always accompanied with at least one space-time 
derivative. 

As discussed, the relevant degrees of freedom are pions (Goldstone bosons) and nucleons. Since the 
interactions of Goldstone bosons must vanish at zero momentum transfer and in the chiral limit (711,,- — )• 0), 
the low-energy expansion of the Lagrangian is arranged in powers of derivatives and pion masses. The hard 
scale is the chiral-symmetry breaking scale, A^^ 1 GeV. Thus, the expansion is in terms of powers of Q/A^ 
where Q is a (small) momentum or pion mass. This is chiral perturbation theory (ChPT). 

The effective Lagrangian can formally be written as, 

^cff = -Ctttt + ^ttN + ^NN + • • • , (1) 

where Ct^tv deals with the dynamics among pions, C-^n describes the interaction between pions and a nucleon, 
and Cnn provides contact interactions between two nucleons. The individual Lagrangians are organized as 
follows: 

^TTTT — ^ir-K + -^iir + ■ • ■ I (2) 
CttN = + ^In + -^i-Af + • • • J (3) 

and 

/-NN — '-NN ' '-NN ' '-NN + • • • 7 [^) 

where the superscript refers to the number of derivatives or pion mass insertions (chiral dimension) and the 
ellipsis stands for terms of higher dimensions. See Ref. [33j for a concise summary of the explicit expressions 
for these Lagrangians. 



3 Nuclear forces from chiral EFT: Overview 
3.1 Chiral perturbation theory and power counting 

As discussed, effective field theories (EFTs) are defined in terms of effective Langrangians which are given by 
an infinite series of terms with increasing number of derivatives and/or nucleon fields, with the dependence 
of each term on the pion field prescribed by the rules of broken chiral symmetry. Applying this Lagrangian 
to a particular process, an unlimited number of Feynman graphs can be generated. Therefore, we need a 
scheme that makes the theory manageable and calculable. This scheme which tells us how to distinguish 
between large (important) and small (unimportant) contributions is chiral perturbation theory (ChPT), and 
determining the power v of the expansion has become known as power counting. 

Nuclear potentials are defined as sets of irreducible graphs up to a given order. The power u oi a, 
few-nucleon diagram involving A nucleons is given in terms of naive dimensional analysis by: 

z/ = -2 + 2A - 2C + 2L + ^ Ai , (5) 

i 

with 

A, = rf, + ^ - 2 , (6) 

where C denotes the number of separately connected pieces and L the number of loops in the diagram; di is 
the number of derivatives or pion-mass insertions and rii the number of nucleon fields (nucleon legs) involved 
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in vertex i; the sum runs over all vertices contained in the diagram under consideration. Note that > 
for all interactions allowed by chiral symmetry. For an irreducible NN diagram ( "two-nucleon force", A = 2, 



In summary, the chief point of the ChPT expansion is that, at a given order there exists only a finite 
number of graphs. This is what makes the theory calculable. The expression {Q / A^Y'^^ provides a rough 
estimate of the relative size of the contributions left out and, thus, of the accuracy at order v. In this sense, 
the theory can be calculated to any desired accuracy and has predictive power. 

3.2 The hierarchy of nuclear forces 

Chiral perturbation theory and power counting imply that nuclear forces emerge as a hierarchy controlled 
by the power v, Fig. [T] 

In lowest order, better known as leading order (LO, z/ = 0), the NN amplitude is made up by two 
momentum-independent contact terms (~ Q°), represented by the four-nucleon-leg graph with a small-dot 
vertex shown in the first row of Fig. [T] and static one-pion exchange (IPE), second diagram in the first 
row of the figure. This is, of course, a rather crude approximation to the two-nucleon force (2NF), but 
accounts already for some important features. The IPE provides the tensor force, necessary to describe the 
deuteron, and it explains NN scattering in peripheral partial waves of very high orbital angular momentum. 
At this order, the two contacts which contribute only in S'-waves provide the short- and intermediate-range 
interaction which is somewhat crude. 

In the next order, v — 1, all contributions vanish due to parity and time-reversal invariance. 

Therefore, the next-to-leading order (NLO) is = 2. Two-pion exchange (2PE) occurs for the first 
time ("leading 2PE") and, thus, the creation of a more sophisticated description of the intermediate-range 
interaction is starting here. Since the loop involved in each pion-diagram implies already z/ = 2 [cf. Eq. ([t])], 
the vertices must have = 0. Therefore, at this order, only the lowest order ttNN and mrNN vertices 
are allowed which is why the leading 2PE is rather weak. Furthermore, there are seven contact terms of 
0{Q^), shown by the four-nucleon-leg graph with a solid square, which contribute in S and P waves. The 
operator structure of these contacts include a spin-orbit term besides central, spin-spin, and tensor terms. 
Thus, essentially all spin-isospin structures necessary to describe the two-nucleon force phenomenologically 
have been generated at this order. The main deficiency at this stage of development is an insufficient 
intermediate-range attraction. 

This problem is finally fixed at order three {ly = 3), next-to-next-to-leading order (NNLO). The 2PE 
involves now the two-derivative ttttNN seagull vertices (proportional to the Ci LECs) denoted by a large solid 
dot in Fig. [l] These vertices represent correlated 2PE as well as intermediate A(1232)-isobar contributions. 
It is well-known from the meson phenomenology of nuclear forces |75l I76| that these two contributions are 
crucial for a realistic and quantitative 2PE model. Consequently, the 2PE now assumes a realistic size 
and describes the intermediate-range attraction of the nuclear force about right. Moreover, first relativistic 
corrections come into play at this order. There are no new contacts. 

The reason why we talk of a hierarchy of nuclear forces is that two- and many-nucleon forces are created 
on an equal footing and emerge in increasing number as we go to higher and higher orders. At NNLO, the 
first set of nonvanishing three-nucleon forces (3NF) occur [3Sll37j, cf. column '3N Force' of Fig. [T] In fact, at 
the previous order, NLO, irreducible 3N graphs appear already, however, it has been shown by Weinberg [TH] 
and others [36l [77l [78] that these diagrams all cancel. Since nonvanishing 3NF contributions happen first at 
order (Q/A^)^, they are very weak as compared to 2NF which start at (Q/A-^)^. 

More 2PE is produced at = 4, next-to-next-to-ncxt-to-leading order (N'^LO), of which we show only a 
few symbolic diagrams in Fig. [T] Two-loop 2PE graphs show up for the first time and so does three-pion 
exchange (3PE) which necessarily involves two loops. 3PE was found to be negligible at this order ^51 124| . 
Most importantly, 15 new contact terms ^ Q"^ arise and are represented by the four-nucleon-leg graph with 
a solid diamond. They include a quadratic spin-orbit term and contribute up to D-waves. Mainly due to 



C = 1), Eq. ([5]) collapses to 




(7) 
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2N Force 



3N Force 4N Force 




NLO 



NNLO 




Figure 1: Hierarchy of nuclear forces in ChPT. Solid lines represent nucleons and dashed lines pions. Small 
dots, large solid dots, solid squares, and solid diamonds denote vertices of index A = 0, 1, 2, and 4, 
respectively. Further explanations are given in the text. 

the increased number of contact terms, a quantitative description of the two-nucleon interaction up to about 
300 MeV lab. energy is possible, at N'^LO (for details, see below). Besides further 3NF, four-nucleon forces 
(4NF) start at this order. Since the leading order 4NF come into existence one order higher than the leading 
3NF, 4NF are weaker than 3NF. Thus, ChPT provides a straightforward explanation for the empirically 
known fact that 2NF » 3NF > 4NF .... 





4 The nucleon-nucleon interaction 

4.1 Definition of the chiral A^A^ potential 

The previous section has provided us with an overview. In this section, we will now discuss the N N 
interaction in more detail. In terms of naive dimensional analysis or "Weinberg counting" , the various 
orders of the irreducible graphs which define the chiral NN potential are given by: 

^Lo = Kf+nl"^ (8) 

^NLO = Vi^o + V^^^ +V^lHv^f (9) 

^^NNLO = V^i^o + Vil^ +V^^^ (10) 

Vn^^o = V^^i^o + V^^ +V^t^ +V^t^ +V^t^ (11) 
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where the superscript denotes the order v of the low-momentum expansion. LO stands for leading order, 
NLO for next-to- leading order, etc.. Contact potentials carry the subscript "ct" and pion-exchange potentials 
can be identified by an obvious subscript. 

The one-pion exchange (IPE) potential reads 

where p ' and p designate the final and initial nucleon momenta in the center-of-mass system and q = p' —p 
is the momentum transfer; Si^2 and ti^2 are the spin and isospin operators of nucleon 1 and 2; gA^ f-K, and 
denote axial-vector coupling constant, the pion decay constant, and the pion mass, respectively. Since 
higher order corrections contribute only to mass and coupling constant renormalizations and since, on shell, 



there are no relativistic corrections, the on-shell IPE has the form Eq. (12 1 up to all orders. 

Multi-pion exchange, which starts at NLO and continues through all higher orders, involves divergent 
loop integrals that need to be regularized. An elegant way to do this is dimensional regularization which 
(besides the main nonpolynomial result) typically generates polynomial terms with coefficients that are, 
in part, infinite or scale dependent 21]. One purpose of the contacts is to absorb all infinities and scale 
dependencies and make sure that the final result is finite and scale independent. This is the renormalization 
of the perturbatively calculated NN amplitude (which, by definition, is the 'WiV potential"). It is very 
similar to what is done in the ChPT calculations of tttt and TriV scattering, namely, a renormalization order 
by order, which is the method of choice for any EFT. Thus, up to this point, the calculation fully meets 
the standards of an EFT and there are no problems. The perturbative NN amplitude can be used to make 
model independent predictions for peripheral partial waves [211 1221 133) . A concise summary of the explicit 
expressions for the 2tt exchange contributions up to order N'^LO can be found in Ref. 



4.2 Regularization and renormalization 

For calculations of the structure of nuclear few and many-body systems, the lower partial waves are the most 
important ones. The fact that in S waves we have large scattering lengths and shallow (quasi) bound states 
indicates that these waves need to be treated nonperturbatively. Following Weinberg's prescription [16] , this 
is accomplished by inserting the potential V into the Lippmann-Schwinger (LS) equation: 

Tip',p) = V{p',p} + f d^p" V{p',p") . T{p",p) , (13) 

J — p" + le 

where Mm denotes the nucleon mass. 

In general, the integral in the LS equation is divergent and needs to be regularized. One way to do this 
is by multiplying V with a regulator function 

V{p\p) ^ V{p',p) e-^P'/^^''" e-(f/^)'" . (14) 

Typical choices for the cutoff parameter A that appears in the regulator are A « 0.5 GeV ^ A^ « 1 GeV. 

It is pretty obvious that results for the T-matrix may depend sensitively on the regulator and its cutoff 
parameter. This is acceptable if one wishes to build models. For example, the meson models of the past [79l 
I75j always depended sensitively on the choices for the cutoff parameters which, in fact, were important for 
the fit of the NN data. However, the EFT approach wishes to be fundamental in nature and not just another 
model. 

In field theories, divergent integrals are not uncommon and methods have been developed for how to deal 
with them. One regulates the integrals and then removes the dependence on the regularization parameters 
(scales, cutoff's) by renormalization. In the end, the theory and its predictions do not depend on cutoffs or 
renormalization scales. 

So-called renormalizable quantum field theories, like QED, have essentially one set of prescriptions that 
takes care of renormalization through all orders. In contrast, EFTs are renormalizcd order by order. 
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As discussed, the renormalization of perturbative EFT calculations is not a problem. The problem is 
nonperturbative renormalization. This problem typically occurs in nuclear EFT because nuclear physics is 
characterized by bound states which are nonperturbative in nature. EFT power counting may be different for 
nonperturbative processes as compared to perturbative ones. Such difference may be caused by the infrared 
enhancement of the reducible diagrams generated in the LS equation. 

Weinberg's implicit assumption [161 113j was that the counterterms introduced to renormalize the per- 
turbatively calculated potential, based upon naive dimensional analysis ("Weinberg counting"), are also 
sufficient to renormalize the nonperturbative resummation of the potential in the LS equation. In 1996, 
Kaplan, Savage, and Wise (KSW) [SIT pointed out that there are problems with the Weinberg scheme if 
the LS equation is renormalized by minimally-subtracted dimensional regularization. This criticism resulted 
in a flurry of publications on the renormalization of the nonperturbative NN problem. See Ref. [8T] for a 
thorough discussion, the bottom line of which is: 

Crucial for an EFT are regulator independence (within the range of validity of the EFT) and a power 
counting scheme that allows for order-by-order improvement with decreasing truncation error. The purpose 
of renormalization is to achieve this regulator independence while maintaining a functional power counting 
scheme. After the comprehensive tries and errors of the past, it appears that there are two renormaliza- 
tion schemes which have the potential to achieve the above goals and, therefore, should be investigated 
systematically in the near future. 

In scheme one, the LO calculation is conducted nonperturbatively (with A — )• oo as in |82| ) and subleading 
orders are added perturbatively in distorted wave Born approximation. Valderrama has started this in S 
waves |83j , but results in higher partial waves are needed to fully assess this approach. Even though at this 
early stage any judgement is speculative, we take the liberty to predict that this approach will be only of 
limited success and utility — for the following reasons. First, it will probably require about twice as many 
counterterms as Weinberg counting and, therefore, will have less predictive power. Second, this scheme may 
converge badly, because the largest portion of the nuclear force, namely, the intermediate-range attraction 
appears at NNLO. Third, as discussed in this force may be problematic (and, therefore, impractical) 
in applications in nuclear few- and many-body systems, because of a pathologically strong tensor force that 
will cause bad convergence of energy and wave functions. Finally, in the work that has been conducted so 
far within this scheme by Valderrama, it is found that only rather soft cutoffs can be used. 

The latter point (namely, soft cutoffs) suggests that one may then as well conduct the calculation nonper- 
turbatively at all orders (up to N'^LO) using Weinberg counting, which is no problem with soft cutoffs. This 
is scheme two that we propose to investigate systematically. In the spirit of Lepage |85j , the cutoff indepen- 
dence should be examined for cutoffs below the hard scale and not beyond. Ranges of cutoff independence 
within the theoretical error are to be identified using 'Lepage plots'. A very systematic investigation of this 
kind does not exist at this time and is therefore needed. However, there is comprehensive circumstantial 
evidence from the numerous chiral NN potentials constructed over the past decade [Ml 1311 [331 [311 [33 [HS] in- 
dicating that this investigation will most likely be a success (cf. also Fig. [2] below) . The potentials discussed 
in the following section are all based upon Weinberg counting. 



4.3 Constructing quantitative chiral A^A^ potentials 
4.3.1 What order? 

As discussed, the NN potential can be calculated up to various orders, cf. Eqs. (|8|-(11), and the accuracy 
increases as the order increases. That triggers the obvious question: To what order do we have to go for an 
accuracy that we would perceive as necessary and sufficient for, e. g., reliable microscopic nuclear structure 
calculations? 

To answer this question, we show in Table[l]the x^/datum for the fit of the world np data below 290 MeV 
for families of np potentials at NLO and NNLO constructed by the Juelich group |89]. The NLO potentials 
produce the very large /datum between 67 and 105, and the NNLO are between 12 and 27. The rate of 
improvement from one order to the other is very encouraging, but the quality of the reproduction of the np 
data at NLO and NNLO is obviously insufficient for reliable predictions. 



8 



Table 1: Columns three and four show the x^/datum for the reproduction of the 1999 np database [57] 
(subdivided into energy intervals) by families of np potentials at NLO and NNLO constructed by the Juelich 
group [M]. The /datum is stated in terms of ranges which result from a variation of the cutoff parameters 
used in the regulator functions. The values of these cutoff parameters in units of MeV are given in parentheses. 
Tiab denotes the kinetic energy of the incident neutron in the laboratory system. 



Jlab bin 


# of np 


— Juelich np potentials — 


(MeV) 


data 


NLO 


NNLO 






(550/700-400/500) 


(600/700-450/500) 


0-100 


1058 


4-5 


1.4-1.9 


100-190 


501 


77-121 


12-32 


190-290 


843 


140-220 


25-69 


0-290 


2402 


67-105 


12-27 



Based upon these facts, it has been pointed out in 2002 by Entem and Machleidt [311 1331 ^^^^ one has 
to proceed to N'^LO. Consequently, the first N'^LO potential was published in 2003 [31]. 

At N^LO, there are a total of 24 contact terms (24 parameters) which contribute to the partial waves with 
L < 2. These 24 LECs are essentially free constants which parametrize the short-ranged phenomenological 
part of the interaction. In Table [2j column 'Q^/N'^LO', we show how these terms are distributed over the 
partial waves. Most important for the improved reproduction of the NN phase shifts (and NN observables) 
at N^LO is the fact that contacts appear for the first time in ZJ-wavcs. _D-waves are not truely peripheral 
and, therefore, IPE plus 2PE alone do not describe them well. The Z3-wave contacts provide the necessary 
short-range corrections to get the ZJ-phases right. Besides this, at N'^LO, another contact is added to each 
P-wave, which leads to substantial improvements, particularly, in ■^Pq and "^Pi above 100 MeV. 

In Table[2j we also show the number of parameters used in the Nijmegen partial wave analysis (PWA93) [90] 
and in the high-precision CD-Bonn potential [88^^. The table reveals that, for S and P waves, the number 
of parameters used in high-precision phenomenology and in EFT at N'^LO are about the same. Thus, the 
EFT approach provides retroactively a justification for the phenomenology used in the 1990's 
to obtain high-precision fits. 

At NLO and NNLO, the number of contact parameters is substantially smaller than for PWA93 and 
CD-Bonn, which explains why these orders are insufficient for a quantitative potential. The 24 parameters 
of N^LO are close to the 30-K used in PWA93 and high precision potentials. Consequently (see following 
sections for details), at N'^LO, a fit of the NN data is possible that is of about the same quality as the one 
by the high-precision NN potentials [SH IM IMIIHH] • Thus, one may perceive N^LO as the order of ChPT 
that is necessary and sufficient for a reliable NN potential. 



4.3.2 A quantitative NN potential at N^LO 



We choose the Idaho N'^LO potential [Mj as example. Note that for an accurate fit of the low-energy pp and 
np data, charge-dependence is important. We include charge-dependence up to next-to- leading order of the 
isospin-violation scheme (NLO, in the notation of Ref. [M])- Thus, we include the pion mass difference in 
IPE and the Coulomb potential in pp scattering, which takes care of the L0 contributions. At order NLO 
we have pion mass difference in the NLO part of TPE, ir^ exchange [SS], and two charge-dependent contact 
interactions of order which make possible an accurate fit of the three different ^5*0 scattering lengths, 
Ojppj Oijini and dnp- 



For the cutoff parameter A of the regulator Eq. (14 1, we choose initially 500 MeV. Within a certain 



reasonable range, results should not depend sensitively on A (cf. discussion in Section 4.2). Therefore, we 
have also made a second fit for A = 600 MeV. 

The fitting procedure starts with the peripheral partial waves because they involve fewer and more 
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Table 2: Number of parameters needed for fitting the np data in phase-shift analysis and by a high-precision A'^A'^ 
potential versus the total number of NN contact terms of EFT based potentials to different orders. 
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fundamental parameters. Partial waves with L > 3 are exclusively determined by IPE and 2PE because 
the N'^LO contacts contribute to L < 2 only. IPE and 2PE at N'^LO depend on the axial-vector coupling 
constant, gA (we use qa = 1-29), the pion decay constant, — 92.4 MeV, and eight low-energy constants 
(LECs) that appear in the dimension-two and dimension-three ixN Lagrangians. The LECs are hsted in 
Table [3) where column WiV pot.' shows the values used for the present N'^LO potential. In the fitting 
process, we varied three of them, namely, C2, C3, and C4. We found that the other LECs are not very effective 
in the NN system and, therefore, we left them at their central values as determined in nN analysis. The most 
influential constant is C3, which — in terms of magnitude — has to be chosen on the low side (slightly more 
than one standard deviation below its t^N determination), otherwise there is too much central attraction. 
Concerning C4, our choice C4 = 5.4 GeV~^ lowers the ^^^2 phase shift (and slightly the ^F-^ bringing it 
into closer agreement with the phase shift analysis — as compared to using the ixN value C4 = 3.4 GeV~^. 
The other F waves and the higher partial waves are essentially unaffected by this variation of C4. Finally, 
the change of ci from its TriV value of 3.28 GeV~^ to 2.80 GeV~^ (our choice) brings about some subtle 
improvements of the fit, but it is not essential. Overall, the fit of all J > 3 waves is very good. 

We turn now to the lower partial waves. Here, the most important fit parameters are the ones associated 
with the 24 contact terms that contribute to the partial waves with L < 2. In addition, we have two 
charge-dependent contacts which are used to fit the three different ^Sq scattering lengths, Cpp, a„„, and a„p. 

In the optimization procedure, we fit first phase shifts, and then we refine the fit by minimizing the 
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Table 3: Low-energy constants applied in the N^LO NN potential (column 'TV A'' pot.'). The q belong to the 
dimension- two ttN Lagrangian and are in units of GeV~^, while the di are associated with the dimension- 
three Lagrangian and are in units of GeV~^. The column 'ttTV analysis' shows values determined from ttN 
data. 







NN pot. 




t:N analysis 


Cl 




-0.81 




-0.81 ± 0.15° 


C2 




2.80 




3.28 ±0.23'' 


C3 




-3.20 




-4.69± 1.34^^ 


_ C4 _ 




5.40 




3.40 ± 0.04^^ 


di + d2 




3.06 




3.06 ±0.21'' 


d_3 




-3.27 




-3.27 ±0.73'' 


d5 




0.45 




0.45 ± 0.42'' 


di4 — di5 




-5.65 




-5.65 ±0.41'' 


"Table 1, Fit 1 of Ref. 


[96]. ''Table 2, Fit 1 of Ref. [97]. 






Table 4: Columns three to five display the x^/datum for the reproduction of the 1999 np database fST] 


(subdivided into energy intervals) by various np potentials. For the chiral potentials, the x^/datum is stated 


in terms of ranges which result from 


a variation of the cutoff parameters used in the rej 


julator functions. The 


values of these cutoff parameters in 


units of MeV are given 


in parentheses. Tiab denotes the kinetic energy 


of the incident nucleon in the laboratory system. 






?iab bin 


# of np 


Idaho 
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Argonne 


(MeV) 


data 


N^LO 133 


N^LO [SB] 


m 






(500-600) 


(600/700-450/500) 




0-100 


1058 


1.0-1.1 


1.0-1.1 


0.95 


100-190 


501 


1.1-1.2 


1.3-1.8 


1.10 


190-290 


843 


1.2-1.4 


2.8-20.0 


1.11 


0-290 


2402 


1.1-1.3 


1.7-7.9 


1.04 



obtained from a direct comparison with the data. We start with pp, since the pp phase shifts and data are 
more accurate than the np ones. The pp fit fixes essentially the 1=1 potential. The I = 1 np potential is 
just the pp one modified by charge-dependence due to nuclcon-mass difference, pion-mass splitting in IPE, 
TTj exchange, and omission of Coulomb. In addition to this, the non-derivative contact in the ^Sq state is 
changed such as to reproduce the np scattering length. The nn potential is the pp one without Coulomb, 
using neutron masses, and fitting the nn scattering length in the ^Sq state with the non-derivative contact. 

The /datum for the fit of the np data below 290 MeV are shown in Table [ij and the corresponding ones 
for pp are given in Table [H] These tables reveal that at N'^LO a x^ /datum comparable to the high-precision 
Argonne Vig |92| potential can, indeed, be achieved. The Idaho N^LO potential [34] with A = 500 MeV 
produces a x^ /datum =1.1 for the world np data below 290 MeV which compares well with the x^ /datum 
= 1.04 by the Argonne potential. In 2005, also the Juelich group produced several N'^LO NN potentials |86l , 
the best of which fits the np data with a x^ /datum = 1.7 and the worse with 7.9 (Table [4|. 

Turning to pp, the x^ for pp data are typically larger than for np because of the higher precision of pp 
data. Thus, the Argonne Vis produces a x^ /datum =1.4 for the world pp data below 290 MeV and the best 
Idaho N'^LO pp potential obtains 1.5. The fit by the best Juelich N'^LO pp potential results in a x^/datum 
= 2.9 and the worst produces 22.3. 

Phase shifts of np scattering from two Idaho (solid and dashed lines) and two Juelich (dash-dotted and 
dotted lines) N^LO np potentials are shown in Figs. ^ The phase shifts confirm what the corresponding x^ 
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Table 5: Same as Table |4] but for pp. 



?iab bin 


# of pp 
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(MeV) 


data 


N^LO 131] 


N^LO [HS] 
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(500-600) 


(600/700-450/500) 




0-100 
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1.0-3.8 
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100-190 
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1.5-1.9 


3.5-11.6 


1.3 


190-290 


851 


1.9-2.7 


4.3-44.4 


1.8 


0-290 


2057 


1.5-2.1 


2.9-22.3 


1.4 



have already revealed. For the low-energy scattering and deuteron parameters, see Ref. |34) . 



5 Few-nucleon forces 

The chiral 2NF discussed in the previous section has been applied in microscopic calculations of nuclear 
structure with, in general, a great deal of success [lTllllllllSllSlinill71llilini[Sni[CTl[51[S3 
However, from high-precision studies conducted in the 1990s, it is well-known that certain few-nucleon 
reactions and nuclear structure issues require 3NFs for their microscopic explanation. Outstanding examples 
are the Ay puzzle of N-d scattering [99 lll00j and the ground state of "'^^B [101 . As noted before, an important 
advantage of the EFT approach to nuclear forces is that it creates two- and many-nucleon forces on an equal 
footing (cf. the overview given in Fig. [T]) . 

For a 3NF, we have ^ = 3 and C — 1 and, thus, Eq. ([5| implies 

ly = 2 + 2L + '^Ai. (15) 

i 

This equation can be used to analyze 3NF contributions order by order. The lowest possible power is 
obviously 1^ — 2 (NLO), which is obtained for no loops (L — 0) and only leading vertices {J^t = 0). This 
3NF happens to vanish [HI |771 [75] . The first non- vanishing 3NF occurs at NNLO. 

The power — 3 (NNLO) is obtained when there are no loops {L = 0) and J^i = 1? i-^-, A.; — 1 for 
one vertex while = for all other vertices. There are three topologies which fulfill this condition, known 
as the two-pion exchange (2PE), IFF, and contact graphs [311137] (Fig. |3|. 

The 2PE 3N-potential is given by 

TA3NF ^ f 9aY I ((7, ■ qi){aj ■ q,j) . 

with qi = Pi' ~ Pi, where pi and p/ are the initial and final momenta of nucleon i, respectively, and 

4ciml 2c3 



f2 f2 

J TT J TT 



|^6"^=r,^a,.[glxg;]. (17) 



There are great similarities between this force and earlier derivations of 2PE 3NFs, notably the 50-year old 
Fujita-Miyazawa [102] . the Tucson-Melbourne (TM) [103j . and the Brazil [104j forces. This is also the 3NF 
advocated by Gerry Brown in Refs. [31 Hj- 
The IFF contribution is 

y?^i = 1^ E #rS(-' • • (18) 
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Figure 2: Neutron-proton phase parameters as described by various chiral potentials at N'^LO. The (red) 
soUd and the dashed curves are calculated from Idaho N'^LO potentials [34] with A = 500 and 600 MeV, re- 
spectively; while the (green) dash-dotted and the dotted curves are based upon Juelich N^LO potentials f5S] 
with cutoff combinations 600/700 and 450/500 MeV, respectively. Partial waves with total angular momen- 
tum J < 2 and laboratory energies up to 300 MeV are displayed. The solid dots and open circles are the 
results from the Nijmegen multi-energy wp phase shift analysis [30] and the VPI/GWU single-energy np 
analysis SM99 [98j, respectively. 



and the 3N contact potential reads 

Vr = E\Y^ryr,. (19) 

The last two 3NF terms involve the two new parameters D and i?, which do not appear in the 2N problem. 
There are many ways to pin these two parameters down. In Ref . |37| , the triton binding energy and the nd 
doublet scattering length ^and were used. One may also choose the binding energies of '^H and ''He [SS] or 
an optimal over-all fit of the properties of light nuclei I66J- Once D and E are fixed, the results for other 
3N, 4N, etc. observables are predictions. 

The 3NF at NNLO has been applied in calculations of few-nucleon reactions [571[5^[571[551[5gi[5gilCTllB^ 
[551 [Ml , structure of light- and medium-mass nuclei [53 [551 [53 EH] , and nuclear and neutron matter [531 IZO] 
with a good deal of success. Yet, the famous ''Ay puzzle' of nucleon-deuteron scattering is not resolved [S71l61| . 
When only 2NFs are applied, the analyzing power in p-^He scattering is even more underpredicted than in 
p-d [40| ,ll05j . However, when the NNLO 3NF is added, the p-'^He Ay substantially improves (more than in -p- 
d) |64j — but a discrepancy remains. Furthermore, the spectra of light nuclei leave room for improvement [55]. 

To summarize, the 3NF at NNLO is a remarkable contribution: It represents the leading many-body force 
within the scheme of ChPT; it includes terms that were advocated by Gerry Brown already some 40 years 
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Figure 3: The three-nucleon force at NNLO. From left to right: 2PE, IPE, and contact diagrams. Notation 
as in Fig. [T] 

ago P|f51; and it produces noticeable improvements in few-nucleon reactions and the structure of light nuclei. 
But unresolved problems remain. Moreover, in the case of the 2NF, we have seen that one has to proceed 
to N'^LO to achieve sufficient accuracy. Therefore, the 3NF at N'^LO is needed for at least two reasons: 
for consistency with the 2NF and to hopefully resolve outstanding problems in microscopic structure and 
reactions. For further discussion of the 3NF (beyond NNLO), see Ref. [81 . 

6 Conclusions and Outlook 

Forty years after Gerry Brown had stressed repeatedly the outstanding importance of chiral symmetry, we 
have finally arrived at a derivation of nuclear forces which takes this symmetry consistently into account. 

The greatest progress occurred during the past 15 years. Key to this development was the realization 
that low-energy QCD is equivalent to a chiral effective field theory which allows for a pcrturbative expansion 
that has become known as chiral perturbation theory (ChPT). In this framework, two- and many-body forces 
emerge on an equal footing and the empirical fact that nuclear many-body forces are substantially weaker 
than the two-nucleon force is explained naturally. 

In this chapter, we have shown how the two-nucleon force is derived from ChPT and demonstrated 
that, at N'^LO, the accuracy can be achieved that is necessary and sufficient for reliable microscopic nuclear 
structure predictions. First calculations applying the N'^LO NN potential [34 in the conventional shell 
model [m 132 US] , the ab initio no-core shell model HH 133 HS] , the coupled cluster formalism [33 |3Hl Ell 
[50l [5TJ [521 153] , ^md the unitary- model-operator approach [54l[55] have produced promising results. 

We also discussed nuclear many-body forces based upon chiral EFT. The 3NF at NNLO has been known 
for a while [3S1 [37] and applied in few-nucleon reactions [371[5ni[SZl[SHl[Ml[Sniinil[Bll[Ml[M], structure of 
light- and medium- mass nuclei [S51 [SHI [13 [SSI , and nuclear and neutron matter [551 [7D] with good success; 
but some open issues remain, which are presently under investigation |81| . 

We are optimistic that the remaining outstanding problems will be resolved within the next few years, 
such that, after 80 years of desperate struggle, we may finally claim that the nuclear force problem is 
essentially under control — last not least, due to a thorough accounting of chiral symmetry. 
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